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A standing mystery in the Standard Model is the unnatural smallness of the strong CP violating 
phase. A massless up quark has long been proposed as one potential solution. A lattice calculation of 
the constants of the chiral Lagrangian essential for the determination of the up quark mass, 2 q 8 — Qs, 
is presented. We find 2as — 05 = 0.29 ± 0.18, which corresponds to mu/md = 0.410 ± 0.036. This 
is the first such calculation using a physical number of dynamical light quarks, Nf = 3. 

PACS numbers: 11.30.Er, 12.39.Fe, 12.38.Gc, 11.15.Ha 


INTRODUCTION 

The nontrivial topological structure of the QCD gauge 
vacuum generates a CP breaking term in the QGD La¬ 
grangian. However, measurements of the neutron electric 
dipole moment have placed a restrictive upper bound on 
this term’s coefficient, 9 < 10“® Q. The unnatural small¬ 
ness of 9 is known as the strong CP problem. 

A massless up quark (m„ = 0) has long been proposed 
as a potential elegant solution to the problem. Chiral 
rotations of the quark mass matrix A4 shift 9, 

9 = 9 + avgdet M. , (1) 

where 0 is a fundamental parameter of the Standard 
Model. However, if = 0, then det Ai — 0 and 
arg det Ai is unphysical, leaving one free to remove the 
CP violating term through a simple field redefinition. 

At leading order, Chiral Perturbation Theory (ChPT) 
appears to rule out the possibility of rriu = 0. The quark 
mass ratios, including mu/md, can be determined using 
ChPT’s LO predictions for the light meson masses. 

At NLO, however, new coefficients appear in the chiral 
expansion which contribute to the meson masses. The 
parameters of ChPT are no longer fully determined by 
experimental data. In fact, it is impossible for ChPT to 
distinguish between the effects of a non-zero up quark 
mass and certain large NLO corrections. This is known 
as the Kaplan-Manohar ambiguity 

Distinguishing between a light and a massless up quark 
requires knowledge of the coefficients of the NLO terms 
in the chiral Lagrangian, the Gasser-Leutwyler (CL) co¬ 
efficients. Specifically, it is the combination of constants 
208 — CTs HI] which appears in Am, the NLO correction 
to the quark mass ratios |^. If this combination falls 
within a certain range, —3.3 < 2a% — 0:5 < —1.5, current 
experimental results can not rule out mu = 0 . 

Various assumptions and phenomenological arguments 
have been used in the past to assemble a somewhat stan¬ 
dard set of values for the coefficients [Q. However, be¬ 
cause these coefficients are physically determined by the 
low-energy non-perturbative behavior of QCD, the lat¬ 
tice offers the best opportunity for a truly first-principles 
calculation. 


PARTIALLY QUENCHED CHIRAL 
PERTURBATION THEORY (pqChPT) 

pqChPT H, ^ is the tool through which one can calcu¬ 
late the CL coefficients on the lattice. pqChPT is distinct 
from standard ChPT in that it is constructed from the 
symmetry of a graded group. This graded group follows 
from the presumed quark content of partially quenched 
QCD (pqQCD): separate valence and sea quark flavors 
in addition to ghost quark flavors, which in perturbation 
theory cancel loop corrections due to valence quarks. 

The Lagrangian of pqChPT up to O(p^) follows, with 
only relevant NLO terms shown. 

C = i-if£i[d^Ud^^U'^]-^—sTr[xU'' + Ux] 

+ U sTr[9^C/a^Ul'] sTr[xC/t -b Ux] 

+ Lg sPY[d^Ud^U^ {xU^ + Ux)] 

- Lq sTr[xt/'^ + Ux] sTr[xf7'l' + Ux] 

- Ls sTr [xU^xU^ + UxUx] + ■■■ , (2) 

where U = exp(2f$//), x = 2/ia“^ diag ({mg, my}), 
contains the pseudo-Goldstone “mesons” of the sponta¬ 
neously broken SU{Nf + Nv\Nv) l'^SU{N j + Nv\Nv )r 
symmetry, and U is an element of that group, ms and 
my refer to the bare lattice quark mass parameters, 
which are related to their dimensionful equivalents via 
the lattice spacing, m!^™ = a~^mx- Three degenerate 
sea quarks were used, Nf = 3, while the number of va¬ 
lence quarks Ny cancels in all expressions, affecting only 
the counting of external states. The constants /, ft, and 
the Li’s are unknown, determined by the low-energy dy¬ 
namics of pqQCD. 

Because the valence and sea quark mass dependence 
of the Lagrangian of pqChPT is explicit and full QCD 
is within the parameter space of pqQCD (my = mg), 
the values obtained for the CL coefficients in a pqQCD 
calculation are the exact values for the coefficients in full 
QCD 0:1- Furthermore, independent variation of va¬ 
lence and sea quark masses allows additional lever arms 
in the determination the the coefficients. Because the 
Nf dependence of the Lagrangian is not explicit, the CL 
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“Starting configuration state: ordered, disordered, or thermal. 
^Denotes a hypercubic blocked ensemble. 

“Spatial volume is 12^ X 16. 



FIG. 1: 16® X 32, P = 5.3, ms = 0.01, hypercubic blocked. 

coefficients are functions of Nf. Thus, it is important to 
use a physical number of light sea quarks, as we have, 
when extracting physical results. 


PREDICTED FORMS 

pqChPT predicts forms for the dependence of the pseu¬ 
doscalar mass and decay constant on the valence quark 
mass, here assuming degenerate sea quarks and degener¬ 
ate valence quarks, and cutting off loops at = Ltt/. 


Ml = (47r/)^zmy|l-I-zmy^2a8 - 05 + ^ y 

-I- — (2mv - ms) In zmy > 
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U = /<i 1 + '^zmv 


zNf 


(mv + ms) In | (my -b ms) | 


( 3 ) 
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where z = 2^a ®/(47r/)^. These forms differ slightly 
from those in pt], as the NLO dependence in the sea 


FIG. 2: 16® X 32, /3 = 5.3, ms = 0.01. The squares and 
dashed curve (diamonds and solid curve) are before (after) 
hypercubic blocking. Data consistent with m^ = 0 would fall 
within the range centered on the dashed-dotted curve. 

quark mass has been absorbed into and /. This is 
allowed as the error due to this change is manifest when 
z appears in the NLO terms, pushing the discrepancy up 
to NNLO. Accounting for these absorbed terms would 
require a systematic study at several sea quark masses. 

The forms above are derived assuming degenerate light 
mesons. However, our use of staggered fermions on some¬ 
what coarse lattices generates significant flavor symmetry 
breaking, and thus a splitting of the light meson masses. 
A detailed analysis of this effect with dynamical quarks 
was recently presented in |]l^ . 

In order to study this error, we applied hypercubic 
blocking to several of our ensembles, using the blocking 
coefficients found in ||^. Because hypercubic blocked 
dynamical quarks were not used when generating these 
ensembles, we are using different Dirac operators for the 
valence and sea quarks. While this procedure may not 
have a clean continuum limit, it is still useful for estimat¬ 
ing the systematic error due to flavor symmetry breaking. 

The ensembles were generated using staggered 
fermions via the inexact HMD R-algorithm |Q, which 
allows one to work at Nf = 3, but involves taking the 
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FIG. 3: 8® X 32, several values of (5 and ms with matched 
lattice spacing, hypercubic blocked. 


I root of the quark determinant. The result is a quark 
action which is non-local at finite lattice spacing, but 
should become local in the continuum limit. 

To determine the value of 2as — 05 , the local pseu¬ 
doscalar correlator was calculated using several valence 
quark masses. The correlators were fit to exponentials, 
while the meson mass and decay constants were simulta¬ 
neously fit to the predicted forms. The results of the fit 
are the values n, f, 05 , and 2as — 05 . Figure displays 
an example of one such fit. 

When fitting, a chiral cutoff point in the valence quark 
mass beyond which one expects pqChPT to break down 
must be chosen. To choose our cutoff, we added my 
values to our fit keeping rav/raK < 1 and keeping 
X^/dof « 1. The cutoff determined for the 16^ x 32 
hypercubic blocked ensemble, in terms of mv/rriK, was 
used for all of the (3 = 5.3, ms = 0.01 ensembles, mx 
denotes an ensemble’s valence quark mass at which the 
mass of the lightest pseudoscalar meson with degenerate 
quarks equals the physical kaon mass. We found 2a8 — 05 
to be very sensitive to our cutoff choice. For the 16^ x 32 
hypercubic blocked ensemble, changing the chiral cutoff 
by ±0.14 in my/m/y shifted 2a8 — <25 by ±0.12. 

Figures |^|| display the quantity 


M^{ms)mv 

M^{mv)ms 


( 5 ) 


suggested by 131. This quantity accentuates the NLO 
terms in M^, as is evident by comparing Figure to Fig¬ 
ures &0. It should be noted, however, that the full forms 
of and /tt, (|) and (|), were used when fitting. When 
calculating Rm, we did not use the simplification shown 
in I p^ , but rather used a full numerator and denomina¬ 
tor. 

For each of the plots, the data points are the result 
of individual fits of the correlator at each valence quark 


FIG. 4: 16® X 32, hypercubic blocked with similar lattice spac- 
ings. The circles and dashed curve (diamonds and solid curve) 
are fully (partially) quenched. The quenched data were fit as 
though partially quenched with rris = 0.01, Nf = 3. 

mass, with jackknife error bars. The curves display the 
result of a simultaneous fit of all the correlators below a 
cutoff in my/mx to the predicted forms of M'^ and /,rj 
and (^), with jackknife error bounds. Solid symbols 
are used below the cutoff, while open symbols are used 
beyond it. Because of our small ensemble sizes, the full 
correlation matrix for many of the ensembles proved to 
be nearly singular. Thus, several of the fits do not fully 
account for data correlation. 

We determined the lattice spacing of our ensembles via 
the static quark potential, using a tree-level corrected 
Coulomb term. The form of this term for hypercubic 
blocked ensembles was taken from [Q. 


RESULTS 

Figure || presents Rm for the 16® x 32 ensemble both 
before and after hypercubic blocking. The application of 
hypercubic blocking altered results significantly, suggest¬ 
ing that the effect of flavor symmetry breaking at these 
lattice spacings is significant. The dashed-dotted curve 
uses the results from the hypercubic blocked ensemble’s 
fit, but replaces the value found for 2 a 8 — 0:5 with one 
consistent with a zero up quark mass. The data clearly 
fall well outside this range. 

To estimate the finite volume error in our result, we re¬ 
peated the calculation in a smaller 12® x 16 x 32 volume, 
holding all other parameters fixed. Fitting this ensem¬ 
ble with the same chiral cutoff as the 16® x 32 ensemble 
resulted in the value 2a8 — 05 = 0.226 ± 0.064. This 
matches our quoted result, suggesting that the finite vol¬ 
ume of our 16® X 32 ensemble is not a large source of 
systematic error. 

Figure presents Rm for several ensembles with a va¬ 
riety of sea quark masses and matched lattice spacing. 
























4 


1 

0 


a -1 
a 

CM 

-2 

-3 


-4 

FIG. 5: Compiled results. The UKQCD Collaboration data 
point is taken from |^, showing their statistical error only. 
This point was calculated using Nf = 2 Wilson fermions, and 
thus its relative agreement suggests small Nf dependence. 
The rightmost point denotes the range allowed by rriu = 0. 

The physical volume of these lattices is similar to the 
physical volume of our 16^ x 32 ensemble. The Columbia 
group has determined several values of the critical /3c and 
rric for the Nf = 3, Nt = 4, finite temperature transition 
I p^ . These ensembles were generated using those bare 
parameters. The trend in 2a^ — 05 with the changing 
sea quark mass is inset in Figure This trend can be 
attributed to the sea quark mass dependence which was 
dropped from Equations ^ and A systematic study 
of this dependence would allow a determination of the 
parameter within the dropped term, 2 a 6 — C 14 [l^ . 

Fully quenched and partially quenched ChPT predict 
different forms for Rm- Thus, one might hope to see the 
effects of quenching though an ensemble’s Rm plot. Fig¬ 
ure H shows the 16^ x 32 partially quenched hypercubic 
blocked ensemble along side a fully quenched ensemble 
with similar lattice spacing. The quenched ensemble was 
analyzed as though partially quenched, using ms = 0 . 01 , 
Nf =3. This procedure does not generate a rigorous 
value for 2 a 8 ~ c« 5 i as this would require the use of fully 
quenched ChPT. However, it could offer insight into the 
magnitude of quenching effects. As Figure ^ shows, the 
effects of quenching are not pronounced. This analysis, 
at its current level of precision, is unable to distinguish 
a fully quenched ensemble from a partially quenched en¬ 
semble of equal lattice spacing. 

The results for 2 q;8 — 05 from each ensemble are com¬ 
piled in Figure While these values vary significantly, 
they do so well outside the range required for a zero up 
quark mass, —3.3 < 2a8 — 0.5 < —1.5. Our quoted re¬ 
sult of 2a8 - 05 = 0.287 ± 0.018"*’"* ± 0.18"^"* comes from 
our hypercubic blocked 16^ x 32 ensemble, where the re¬ 
ported systematic error is the result of adding in quadra¬ 
ture the determined effects of shifting the chiral cutoff 
(±0.12), hypercubic blocking (±0.05), doubling the lat¬ 
tice spacing (± 0 . 11 ), and reducing the lattice volume 
(±0.06). Assuming Dashen’s rule Q, this corresponds 
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to Am = -0.0897 ±0.0313 and mu/ma = 0.484 ±0.027, 
where the quoted error arises primarily from the sys¬ 
tematic error of our measurement. The error from ex¬ 
perimental input is negligible and the size of NNLO 
corrections to Am are assumed to be on the order of 
A^. Using the values for the electromagnetic contri¬ 
butions to the light meson masses from |j^ in place of 
Dashen’s rule results in Am = —0.0898 ± 0.0313 and 
iTiu/md = 0.410 ± 0.036. This can be compared to 
previous calculations in the literature, which have given 
iTiu/md = 0.553 ± 0.043 and mu/md = 0.46 ± 0.09 
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